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不 同 地 理 种 群 遗传 分 化 及 基因 流 研 


1 1 1 ,x 1 2 
李 善 , 张 闫 , 王 振 营 '", 何 康 来 , 王 强 
(1. 中 国 农业 科学 院 植物 保护 研究 所 , 植物 病虫害 生物 学 国家 重点 实验 室 , 北京 100193 ; 
2. 北京 城市 学 院 生 物 技术 学 部 , 北京 100083 ) 


摘要 : 本 文通 过 对 中 国 亚 洲 玉 米 贱 Osirinia furnacalis (Guenée)17 个 地 理 种 群 的 线粒体 CO 下 基因 进行 序列 分 析 ， 对 

中 国 亚洲 玉米 蜡 种 群 间 的 遗传 分 化 程度 和 基因 流水 平 进 行 了 初步 研究 。 在 获得 的 413 条 序列 样本 中 共 发 现 了 34 个 
变异 位 点 、35 种 单 倍 型 。 总 体 单 倍 型 多 样 性 指数 Hd 为 0. 811, 种 群 内 单 倍 型 多 样 度 在 0.424 ~ 0. 862 范围 内 。17 个 

种 群 间 的 平均 基因 流 (Nm) 为 2.02。 总 群体 的 固定 系数 Fsi 为 0.234。AMOVA 分 子 方差 分 析 结 果 表 明 中 国 玉 米 蜡 的 

遗传 分 化 主要 来 自 种 群 内 部 (76.45% ) 。 各 种 群 的 Tajima's D 值 中 性 检验 符合 中 性 突变 , 说 明 中 国 亚 洲 玉 米 赋 在 历 

史上 没有 出 现 群 体 扩张 , 群体 大 小 稳定 。 各 地 理 种 群 的 遗传 距离 与 地 理 距 离间 不 具有 显著 的 相关 性 。 各 地 理 种 群 中 

的 单 倍 型 在 系统 发 生 树 上 散布 在 不 同 的 分 布 群 中 , 缺乏 明显 的 地 理 分 布 格局 。 
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Genetic differentiation and gene flow among different geographical 

populations of the Asian corn borer, Ostrinia furnacalis ( Guenee ) 

(Lepidoptera: Crambidae ) in China estimated by mitochondrial CO I 


gene sequences 
LI Jing , ZHANG Ying , WANG Zhen-Ying ’*, HE Kang-Lai , WANG Qiang” (1. State Key Laboratory 
for the Biology of the Plant Diseases and Insect Pests, Institute of Plant Protection, Chinese Academy of 
Agricultural Sciences, Beijing 100193, China; 2. College of Biotechnology, Beijing City University, 
Beijing 100083 ，China ) 

Abstract: To investigate the genetic differentiation and gene flow among different populations of the Asian 
corn borer, Ostrinia furnacalis (Guenée), in China, we analyzed CO I[ sequences of mitochondrial DNA in 
17 geographical populations. In the total 413 CO IL gene sequences obtained, 34 variable sites and 35 
haplotypes were observed. The mean haplotype diversity of all populations was 0. 811, whereas it varied 
from 0. 424 to 0. 862 within population. The average gene flow (Nm) among the 17 populations was 2.02. 
The overall Fst value was 0.234. Molecular variance (AMOVA ) analysis indicated that a high proportion of 
the total genetic variance was attributable to variations within populations (76.45% ). The Tajima’s D value 
neutrality test results suggest that for the Asian corn borer, there might be no population expansion, so their 
population sizes were steady. The genetic distance was not significantly correlated with the geographical 
distance between populations. The phylogenetic tree revealed no obvious geographical pattern for the 
prevalence and distribution of the haplotypes. 
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亚洲 玉米 坚 Ostrinia furnacalis ( Guenée ) 在 系统 秆 野 蜡 属 Ostrinia, 是 世界 性 的 玉米 害虫 , 给 玉米 生 
分 类 上 属 鲜 怒 目 (Lepidoptera) 草 蝶 科 (Crambidae) 产 造成 严重 的 经 济 损失 。 其 寄主 范围 广 , 除 主要 为 
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害 玉米 外 , 还 为 害 高 梁 、 谷 子 、 棉 花 和 甘 基 等 多 种 重 
要 的 粮食 作物 和 经 济 作物 ( 王 振 膏 等 , 2000)。 亚 洲 
玉米 螟 在 我 国 分 布 十 分 广泛 , 其 分 布 区 跨越 了 很 广 
的 经 纬度 地 理 范 围 , 各 地 的 目 然 环 境 、 植 被 类 和 群 、 耕 
作 制 度 、 寄 主 植物 和 种 植 的 玉米 品种 差异 较 大 。 玉 
米 蜡 在 不 同 的 地 理 分 布 区 具有 不 同 的 化 性 , 我 国 自 
北向 南 1 年 可 发 生 1 ~7 代 ( 胡 明 峻 和 孙 伯 欣 ， 
1979) 。 亚 洲 玉 米 晓 成 虫 扩 散 能 力 强 , 且 不 排除 部 
分 个 体 随 气流 远 距 离 迁 飞 的 可 能 ( 王 振 营 等 ， 
1994) 。 由 于 玉米 螟 成 忠 具有 相对 较 强 的 飞行 扩散 
能 力 , 这 种 非 远 距离 的 扩散 知 在 某 种 条 件 下 连续 或 
渐次 地 发 生 , 经 过 奉 干 代 ( 年 ) 后 , 则 有 可 能 使 得 相 
邻 的 \ 不 同步 的 生态 区 域 的 种 群 之 间 发 生 基 因 交 流 。 

对 于 大 多 数 物 种 , 由 于 环境 因子 的 不 对 称 性 ， 
以 及 生物 行为 习性 的 差异 , 可 能 导致 基因 流 的 不 随 
机 性 , 使 区 域 种 群 间 出 现 不 同 程度 的 遗传 差异 。 种 
群 遗传 结构 多 样 性 的 产生 , 为 种 群 间 进一步 分 化 其 
至 形成 新 的 物种 提供 了 机 会 (Avise，1992)。 随 着 
分 子 系统 地 理学 (molecular phylogeography ) 的 兴起 
与 发 展 , 使 得 以 分 子 生 物 学 技术 为 手段 、 在 分 子 水 
平 上 研究 物种 间 及 物种 内 不 同 种 群 形成 现 有 分 布 格 
局 的 历史 原因 和 演化 过 程 成 为 可 能 (Avise，1998 ; 
王 静 等 , 2001) 。 由 于 动物 线粒体 DNA (mtDNA ) 具 
有 母系 遗传 、 进 化 速率 快 , 拷贝 数 多 的 特点 , 和 常 作 
为 研究 物种 间 系 统 进 化 及 种 下 微 进 化 的 有 效 遗 传 标 
记 (Arvise et al., 1987; Harrison，1989 ; Arvise， 
1994) 。 在 昆虫 系统 地 理学 研究 中 , 细胞 色素 C 氧 
化 酶 亚 基 CO I 、CO 本 基因 , 以 及 NADH 氧化 还 原 
酶 亚 基 ND1、NDS 基因 等 常 作为 分 析 亲 绿 关 系 较 近 
的 种 、 亚 种 及 地 理 种 群 之 间 系 统 关 系 的 分 子 工具 
( 王 备 新 和 杨 寻 基 , 2002; 卜 云 和 郑 哲 民 , 2005 ) 。 

本 研究 利用 线粒体 CO I 基因 作为 分 子 标记 ， 
对 中 国 亚洲 玉米 螟 17 个 地 理 种 群 的 遗传 结构 及 基 
因 流 水 平 进行 了 分 析 , 探讨 了 亚洲 玉米 坚 各 地 理 种 
群 之 间 的 遗传 分 化 程度 , 并 初步 推测 了 其 种 群 演化 
格局 可 能 的 形成 原因 ,为 今后 针对 中 国 不 同 亚 洲 玉 
米 虹 发 生 区 制定 有 效 种 群 控 制 策略 及 转基因 抗 虫 玉 
米 的 抗 性 治理 策略 提供 科学 依据 。 同 时 , 本 研究 结 
果 还 将 为 亚洲 玉米 蝶 是 否 会 进行 远 趾 离 迁 飞 的 问题 
在 分 子 水 平 上 提供 一 个 佐证 。 


1 材料 与 方法 


1.1 供 试 虫 源 
本 试验 中 所 用 的 亚洲 玉米 螟 种群, 分别 于 2007 


年 10 月 -2009 年 9 月 采集 目 全 国 17 个 不 同 地 点 的 
玉米 田 。 种 群 样本 、 采 集 地 点 分 布 详情 如 表 1 所 示 。 
采 回 后 单 头 幼虫 置 于 离心 管 中 ，- 80Y 保存 至 DNA 
提取 。 
1.2 基因 组 DNA 提取 

单 头 玉 米 蝶 幼虫 在 1 x TE (Tris-EDTA, pH = 
8.0) 缓 冲 液 中 冰 浴 匀 浆 , 应 用 组 织 / 细 胞 基因 组 
DNA 快速 提取 试剂 盒 (Biomed ) 提取 基因 组 DNA 
(具体 提取 过 程 按照 试剂 盒 说 明 书 步骤 进行 操作 ) 。 
提取 的 DNA 样品 取 3 ~4 pL 通过 1.0% 琢 脂 糖 北 
胶 电 泳 来 检测 提取 质量 , 检测 合格 后 保存 于 -207 
备用 。 
1.3 COEI 基 因 了 PCR 扩 增 及 序列 测定 

根据 GenBank 中 登录 的 亚洲 玉米 坚 线 粒 体 基 
因 组 序列 信息 (GenBank 登录 号 : AF467260) , 通过 
软件 Primer Premier 5.0 设计 用 于 扩 增 CO 本 基因 全 
长 特异 性 引物 ( 扩 增 产物 包括 CO 基因 两 端的 侧 
轰 序 列 )。 引 物 序 列 : COI_F ( 正 问 引物 ): 5 
CCACCGGCAGAACATTCATAT-3’; COI_R (反问 
引物 ): 5'’-GACCATTACTTGCTTTCAGTCATC-3', 扩 
增 目 标 片 段 为 824 bp。 

PCR 反应 总 体系 为 25 pL, 含 2.5 pL 10 x Ex 
Taq Buffer (Mg plus) , 0.1 mmol/L dNTPs, 正 / 反 
癌 引 物 各 0.5 pmol/L, 1.0 pL 模板 DNA, 0.5U Ex 
Taq” DNA 聚合 酶 (大 连 宝 生物 TaKaRa) 。PCR 扩 
增 程序 : 94Y 预 变性 3 min, 94%C 30 s,，60% 退火 1 
min, 72% 延伸 1 min, 反应 35 个 循环 , 72% 延伸 10 
min。 反 应 结束 后 , 取 3 pL PCR 产物 于 1.5% 的 琉 

§ 糖 凝 胶 电泳 检测 。 电 泳 结 果 在 Dolphin-DOC 

( Wealtec Corp. ) 凝 胶 成 像 系 统 下 进行 观察 。 

电泳 检测 确认 成 功 扩 增 后 , 将 PCR 产物 经 纯 
化 后 在 ABI 3730 XL 型 自动 分 析 仪 上 进行 测序 ， 
测序 实行 正 反 双 问 测序 , 以 提高 所 得 序列 的 精 
确 度 。 
1.4 序列 整理 及 数据 分 析 

测序 结 采 通过 Chromas 软件 读 取 , 通过 序列 整 
理 去 除 CO 下 基 因 两 端 侧 履 序列。 将 所 得 序列 与 
GenBank 中 登陆 的 亚洲 玉米 蜡 线 粒 体 CO 本 基因 
(GenBank 登录 号 : FJ435440) 进行 序列 比 对 , 同时 
在 NCBI 网 站 中 进行 BLAST (http://blast. ncbi. 
nlm. nih. gov/Blast. cgi) ， 以 明确 得 到 的 序列 确 为 目 
标 基 因 。 利 用 Clustal X 1.8 (Thompson et al., 1997) 
进行 多 序列 同 源 比 对 , 并 辅 以 人 工 校对 。 
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表 1 亚洲 玉米 晶 地 理 种 群 采集 信息 及 供 试 个 体 数 量 
Table 1 Collecting data and number of assayed individuals of Ostrinia furnacalis in this study 
序号 种 群 代码 采集 地 点 地 理 坐 标 采集 时 间 ( 月 /年 ) 分 析 个 体 数 
No. Population Sampling location Geo-coordinates Collection date Number of 
code ( month/year) individuals assayed 
1 KSh 龙 江 现 48.21°N, 125.66°FE 10/2008 24 
Keshan, Heilongjiang 
十 林 从 主 崔 
2 GZL 所 林 公主 IR 43.31°N, 124.49°F 10/2008 24 
Gongzhuling, Jilin 
下 
3 ChF 站 村 三 赤 几 42.17°N, 118.5S8°FE 10/2008 24 
Chifeng, Inner Mongolia 
山 4 J 中 
4 XZh 图 做 放 38.24°N, 112.43°E 09/2008 24 
Xinzhou, Shanxi 
河北 衡 z 
3 HSh 9 北 重水 37.44°N, 115.42°E 10/2007 24 
Hengshui, Hebei 
山 ~ A ) 中 
6 JZh 东 胶 州 36.17°N, 120. 00°E 10/2007 24 
Jiaozhou, Shandong 
7 GG 甘 址 甘谷 34.84°N, 105.32°FE 10/2008 24 
Gangu, Gansu 
8 MX 陕西 由 县 34.41°N, 107.09°E 09/2009 24 
Meixian, Shaanxi 
河 丙 洛阳 
9 LoY 南 活 34.41°N, 112.27°E 10/2008 24 
Luoyang, Henan 
Lr ay J | 
10 SuZ 妇 知 得 33.38°N, 116. 58°E 09/2008 27 
Suzhou, Anhui 
江苏 J 中 
11 YZh 苏 场 州 32.23°N, 119.26°E 11/2008 24 
Yangzhou, Jiangsu 
新 疆 昌 二 
12 ChjJ 新 疆 虽 于 44.02°N, 87.18°E 10/2007 24 
Changji, Xinjiang 
四 川 郭 县 
13 PX 30.79°N, 103. 87°E 10/2008 27 
Pixian , Sichuan 
浙江 东阳 
14 DY 29.16°N, 120. 14°E 11/2008 24 
Dongyang, Zhejiang 
云南 昆 
15 KM 南 昆明 25.13°N, 102.76°E 08/2009 24 
Kunming, Yunnan 
局 J | 
16 GZh 三 东 4 23.08°N, 113.14°E 11/2008 24 
Guangzhou, Guangdong 
海南 三 亚 
17 SYa 海南 18.14°N, 109.31°E 02/2008 27 


Sanya, Hainan 


应 用 软件 DnaSP 5. 0 (Librado and Rozas, 
2009 ) 计算 种 群 内 及 种 群 间 的 变异 百分比 、 单 伴 型 
多 样 性 友 、 核 苷 酸 多 样 性 指数 +、 种 群 间 遗传 分 化 
程度 参数 Fst 和 Cs 久 种 群 间 核 昔 酸 睹 异 度 Pxy、 核 苷 
酸 平均 差异 数 KE (Jukes-Cantor) 等 分 子 遗 传 学 参数 ， 
计算 种 群 间 的 基因 流 Nm， 并 进行 中 性 检验 
( Kimura ,1983 ; Tajima，1983 , 1989 ) 。 根 据 经 纬度 


计算 出 各 采集 地 两 两 之 间 的 地 理 距 离 。 通 过 MEGA 
4.0 软件 (Tamura et al.，2007 ) 基于 Kimura 2- 
Parameter ( K2P) 模型 ,以 家 春 Bombyx mor CO 工 
基因 作为 外 群 (AB083339 ) ， 采 用 NJ 邻接 法 
(Neighbor-joining) 构建 单 倍 型 系统 发 生 树 ( Tamura 
et al., 2004)。 系 统 树 各 分 支 的 置信 和 度 (bootstrap) 
均 进 行 1 000 次 重复 检验 。 利 用 Arlequin 3. 11 软件 
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(Excoffier et oj/.，2005 ) 计算 种 群 间 的 遗传 距离 ,并 
进行 分 子 变异 分 析 (AMOVA ) 。 为 了 检验 种 群 间 的 
遗传 分 化 是 否 与 种 群 间 的 地 理 距离 存在 相关 关系 ， 
将 成 对 种 群 间 的 遗传 距离 [Fsi/(1 - Fst)] 和 矩阵 与 
各 种 群 采 样 点 间 的 地 理 距离 (km) 的 自然 对 数 用 ZT 
软件 包 检 验 两 矩阵 的 相关 性 (Bonnet and Van de 
Peer, 2002 ) 。 


2 结果 与 分 析 


2.1 亚洲 玉米 旺 线 粒 体 CO II 基因 核 苷 酸 多 态 性 
及 单 们 型 分 析 

亚洲 玉米 蜡 线 粒 体 CO 了 [基因 全 长 682 bp, 在 
17 个 种 群 共 413 条 CO 本 基因 全 长 序列 中 共 发 现 变 
异 位 点 34 个 , 约 占 分 析 位 点 总 数 的 4.99% , 其 中 
单一 多 态 位 点 11 个 , 简约 信息 位 点 23 个 。 在 所 测 
序列 中 没有 发 现 碱 基 缺 失 或 插入 现象 。CO 工 基 因 
G +C 平均 含量 为 22.8% ,说 明 该 基因 富 含 A +T， 
表现 出 明显 的 碱 基 偏 从 性。 据 报 道 , 昆虫 线粒体 
CO 荆 基 因 在 不 同 的 目 中 A +T 含 量 范 围 不 同 ， 鳞 翅 
目 等 完全 变态 类 昆虫 中 A +T 含 量 最 高 , 达 75% 左 
右 , 与 本 研究 所 得 结果 相符 (Jermiin and Crozier， 
1994; Frati et al., 1997 ) 。 在 所 有 的 核 苷 酸 蔡 换 中 ， 
转换 (transition ) 占 96.68% , 其 中 A 与 G 之 间 的 转 
换 占 45.03% ，C 与 T 之 间 的 转换 占 51.65% ; 颠 换 
(transversion ) 发 生 频 率 很 低 , 仪 占 3.32% ; 总 体 转 
换 / 颠 换 偶 倚 率 丸 值 为 13.31。 据 研究 发 现 , 转换 / 
颠 换 率 在 不 同 进化 中 离 的 分 类 阶 元 之 间 具 有 明显 差 
异 , 通常 亲缘 关系 较 近 的 分 类 阶 元 之 间 核 苷 酸 蔡 换 
表现 为 转换 率 大 于 丰 换 率 , 而 在 较为 远 绿 的 分 类 阶 
元 间 核 甘酸 蔡 换 则 以 颠 换 为 主 (Simon ei al., 1994; 
Frati et wj.，1997) 。 本 研究 发 现 亚洲 玉米 蜡 种 内 线 
粒 体 CO 开 基 因 核 苷 酸 多 态 位 点 变异 类 型 以 转换 为 
主 , 验证 了 上 述 规律 。 

实验 中 获得 的 413 条 CO 本 基因 全 长 序列 共 定 
义 了 35 种 单 倍 型 ， Hapl ~ Hap35 (GenBank 登录 
号 : HM368402 ~ HM368436 ) , 其 中 Hap3 为 共享 单 
倍 型 , 占 所 有 样本 的 37.53% , 其 他 在 种 群 中 出 现 
频率 较 高 的 单 倍 型 依次 为 Hap2 (14.04% )、Hap8 
(9.20% ) 、Hap6 (8. 47% )、Hapl$ (8. 23% ) 、Hapl9 
(7.75% )。 为 有 2 种 单 倍 型 均 为 相应 种 群 所 特有 。 
2.2 COI 基 因 各 单 倍 型 间 的 系统 进化 分 析 

采用 邻接 法 (NJ) 构建 亚洲 玉米 虹 CO 本 基 因 
35 种 单 倍 型 之 间 的 系统 树 , 网 中 仅 显 示 出 > 30% 


的 bootstrap 置信 和 度 值 (图 1)。 由 系统 树 可 以 看 出 ， 
单 倍 型 Hap19，Hap32 和 Hap34 单独 聚 为 1 文 
(bootstrap 置信 度 为 99% ) 。 其 中 单 倍 型 Hap19 出 
现 于 甘谷 、. 郭 县 昆明、 宿州 忻州、 扬州 和 三 亚 7 个 
种 群 内 ,而 聚 类 更 紧密 的 Hap32 和 Hap34 仅 出 现在 
忻州 种 群 中 。 除 这 3 种 单 倍 型 外 ,其 余 分 别 有 18 
种 和 14 种 单 倍 型 聚 为 2 大 分 支 (bootstrap 置信 和 度 
为 79% ) 。 
2.3 亚洲 玉米 旦 不 同 地 理 种 群 的 遗传 多 样 性 .种 群 
间 的 基因 流 及 遗传 分 化 分 析 

对 亚洲 玉米 虹 17 个 地 理 种 群 COE 基 因 序 列 采 
用 Tajima's D 值 进 行 中 性 检验 ( 表 2) , 检测 结果 均 


Hap3 
Hap18 
Hap33 
Hap7 
Hap26 
Hapll 
Hap25 
Hap9 
Hap28 
Hap35 
Hap12 
Hap4 
60 Hap31 
Hap20 
Hapl4 
Hapl3 
39 Hap27 
Haps 
79 Hapl0 
45 Hap%8 
42 Hap17 
Hap21 
Hapl6 
Hap2 
44 Hapls 
Hap22 
Hap29 
Hap23 
Hap30 
了 ap6 
31 Hapl 
43 Hap24 
Hap19 
99 Hap32 
79 Hap34 


Bombyx mori 
图 1 利用 邻接 法 构建 的 亚洲 玉米 蜡 不 同 地 理 
种 群 CO 开 单 倍 型 间 的 系统 树 
Neighbor joining phylogenetic tree of haplotypes 
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Fig. 1 
of different geographic populations of Ostrinia furnacalis 


based on CO gene sequence 
以 家 奋 为 外 群 。Bombyx mori was used as the outgroup. 
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不 显著 , 除 GG (甘谷 ) 和 LoY (洛阳 )2 个 种 群 0. 10 

>P>0.05 以 处, 其余 种 群 均 为 已 >0.10, 说 明 中 
国 的 玉米 量 种 群 在 过 去 没有 出 现 群体 扩张 和 持续 增 
长 模式 , 群体 大 小 保持 相对 稳定 状态 。 

各 种 群 CO 卫 单 倍 型 多 样 性 、 核 车 酸 多 样 性 及 
Tajima’s D 中 性 检验 结果 如 表 2 所 示 。 由 结果 可 以 
看 出 , 所 测 的 玉米 曝 17 个 地 理 种 群 的 总 体 单 倍 型 
多 样 性 指数 Ha 为 0. 811, 种 群 间 核 苷 酸 多 样 性 指 


数 7 为 0.00659。 种 群 内 单 倍 型 多 样 度 为 0. 424 ~ 
0.862, 平均 为 0.712, 其 中 单 倍 型 多 样 度 最 高 的 种 
群 为 KM (昆明 ) 种 群 , 而 YZh (扬州 ) 种 群 在 所 检 
测 的 种 群 中 单 倍 型 多 样 度 最 低 。 种 群 内 核 甘酸 多 样 
度 范 围 在 0. 00308 ~ 0. 01015 之 间 , 平均 为 
0.00509, 其 中 核 背 酸 多 样 度 最 高 和 最 低 的 种 群 分 
别 为 GG (甘谷 ) 和 DY (东阳 ) 种 群 。 


表 2 亚洲 玉米 蜡 不 同 地 理 种 群 CO 基因 单 倍 型 、 核 彰 酸 多 样 性 分 析 及 Tajima’”s D 中 性 检验 
Table 2 COI haplotype diversity, nucleotide diversity and Tajima’s D neutrality test 


in different geographic populations of Ostrinia furnacalis 





种 群 代码 单 倍 型 多 样 性 核 苷 酸 多 样 性 核 背 酸 平均 差异 数 中 性 检验 显著 性 
Population Haplotype diversity Nucleotide diversity Average number of Tajima’sD 

code (Hd) (7) nucleotide difference (hk) Statistial significanee 
KSh 0.630 0. 00403 2.750 1.4597 Not significant, 已 >0. 10 
GZL 0.731 0. 00388 2.648 1.2528 Not significant, P >0. 10 
ChF 0. 830 0.00473 3.228 1.1127 Not significant, 已 >0. 10 
XZh 0. 500 0. 00636 4.341 0. 0465 Not significant, 已 >0. 10 
HSh 0.779 0. 00436 2.972 0.7254 Not significant, 已 >0. 10 
JZh 0.779 0. 00369 2.514 0. 5574 Not significant, 已 >0. 10 
CC 0.714 0.01015 6.920 1.8649 Not significant, 0.10>P>0.05 
MX 0.708 0.00413 2.814 1. 5296 Not significant, P >0. 10 
LoY 0.754 0. 00431 2.938 1.7739 Not significant, 0.10>P>0.05 
SuZz 0.709 0. 00526 3.590 _0.4693 Not significant, 已 >0. 10 
YZh 0. 424 0. 00762 5.196 0.7571 Not significant, 已 >0. 10 
CHJ 0.736 0. 00325 2.217 _0.2621 Not significant, 已 >0. 10 
PX 0.795 0. 00563 3.838 _0.7838 Not significant, P >0. 10 
DY 0.632 0. 00308 2. 103 _0.4562 Not significant, P >0. 10 
KM 0. 862 0. 00694 4.732 0. 1414 Not significant, 已 >0. 10 
GZh 0. 808 0. 00379 2. 583 1.1819 Not significant, P >0. 10 
SYa 0.709 0. 00526 3. 590 _0.4693 Not significant, 已 >0. 10 

总 群体 Total 0.811 0. 00659 4.494 _0.4134 Not significant, P >0.10 


Fst( -statistic ) 可 表示 不 同 群 体 间 等 位 基因 频 
率 的 变异 , 是 反映 群体 进化 历史 的 重要 参数 ,可 在 
一 定 程 度 上 揭示 种 群 间 基 因 流 和 遗传 漂 变 的 程度 。 
而 基因 流 则 可 以 揭示 出 群体 间 可 能 的 基因 渗透 及 影 
响 遗 传 分 化 的 遗传 现象 。DnaSP 5.0 软件 计算 结 
显示 , 17 个 种 群 间 的 基因 流 平均 为 2.02, 总 体 遗 传 
分 化 系数 Gst 为 0. 110, 总 群体 的 固定 系数 Fst 为 





0.235。 种 群 间 的 核 苷 酸 平均 差异 数 Kxy 在 2. 105 

~10. 894 之 间 , 平均 为 4.556; 核 背 酸 攻 异 度 Dxy 
在 0.0031 ~0.0160 之 间 , 平均 为 0.0067( 表 3) ; 种 
群 间 Fst 值 的 范围 为 -0.0398 ~0.6420, 其 中 YZh 
(扬州 ) 种 群 与 其 他 种 群 的 遗传 分 化 程度 最 
高 ( 表 4)。 
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表 3 亚洲 玉米 量 成 对 种 群 间 的 核 背 酸 平 均 差异 数 (Kxy) (上 三 角 ) 与 核 苷 酸 收 异 度 (Dzxy) (下 三 角 ) 
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Table 3 Pairwise average number of nucleotide differences (Kxy) (above the diagonal) and nucleotide divergence (Dxy) 


KSh 
GZL 
ChF 
XZh 
HSh 
JZh 


LoY 
SuZ 
YAh 
chJ 
PX 
DY 


GZh 
SYa 


KSh 
GZL 
ChF 


HSh 


8 


SUZ 
YAh 
chJ 
PX 
DY 


GZh 
SYa 


KSh 


0.0044 
0.0043 
0.0071 
0.0044 
0.0043 
0.0097 
0.0043 
0.0045 
0.0048 
0.0159 
0.0045 
0.0051 
0.0044 
0.0056 
0.0042 
0.0048 


表 4 


KSh 


6.31 
6.75 
7.34 
7.29 
7.26 
7.72 
7.69 
7.55 
7.31 
7.33 
8.01 
7.90 
7.68 
8. 10 
8.01 
8.21 


GZL 
2.987 


0.0044 
0.0076 
0.0045 
0.0046 
0.0087 
0.0044 
0.0044 
0.0047 
0.0152 
0.0048 
0.0048 
0.0047 
0.0055 
0.0042 
0.0047 


( below the diagonal) between different geographic populations of Ostrinia furnacalis 


ChF 
2.941 
2.980 


0.0052 
0.0047 
0.0046 
0.0095 
0. 0045 
0.0046 
0.0050 
0.0149 
0.0048 
0.0052 
0.0047 
0.0058 
0.0044 
0.0050 


XZh 
4.813 
3.174 
5.063 


0.0065 
0.0056 
0.0105 
0.0062 
0.0064 
0.0077 
0.0133 
0.0052 
0.0073 
0.0052 
0. 0080 
0.0078 
0.0077 


HSh 
3.031 
3.085 
3.176 
4.406 


0. 0040 
0. 0096 
0.0041 
0.0042 
0.0049 
0.0156 
0. 0040 
0. 0050 
0. 0039 
0.0057 
0.0046 
0.0049 


JZh 
2.934 
3.105 
3.135 
3.854 
2.732 


0. 0095 
0. 0038 
0. 0040 
0. 0050 
0.0152 
0.0034 
0.0049 
0. 0033 
0. 0097 
0. 0047 
0. 0050 


GG 
6. 580 
3.920 
6.469 
7.146 
6. 540 
6.472 


0.0093 
0.0094 
0.0095 
0.0116 
0.0095 
0.0092 
0.0095 
0.0095 
0.0094 
0.0095 


MX 
2.944 
2.9/4 
3.085 
4.232 
2.788 
2.612 
6.359 


0.0041 
0.0048 
0.0153 
0.0038 
0.0049 
0.0037 
0.0056 
0.0045 
0.0048 


LoY 
3.045 
3.031 
3. 160 
4.389 
2. 850 
2.722 
6.417 
2.766 


SuZ 
3.295 
3.209 
3. 390 
5.241 
3.309 
3.383 
6.469 
3.269 
3.289 
0.0048 
0.0154 
0.0040 
0.0049 


0.0155 
0.0051 
0.0054 
0.0039 0.0050 
0.0057 0.0060 
0.0046 0.0048 
0.0048 0.0051 


YZh 
10. 854 
10. 384 
10.785 
9.076 
10.667 
10. 358 
7. 896 
10. 406 
10. 500 
10. 540 


0.0149 
0.0148 
0.0149 
0.0146 
0.0160 
0.0155 


ChJ 
3.056 
3.254 
3.260 
3.556 
2.717 
2.319 
6. 503 
2.594 
2.726 
3.506 
10.177 


0.0049 
0.0031 
0.0058 
0. 0050 
0.0051 


PX 
3.478 
3.245 
3.528 
4.966 
3.419 
3.326 
6.279 
3.304 
3.372 
3.711 
10. 113 
3.301 


0.0049 
0.0061 
0.0051 
0.0048 


DY 
2.980 
3.219 
3. 190 
3.514 
2.641 
2.268 
6.507 
2.322 
2.039 
3.401 
10. 192 
2. 109 
3.316 


KM 
3.854 
3.721 
3.944 
5.444 
3.913 
3.878 
6. 490 
3.819 
3.875 
4.090 
9.990 
3.951 
4.170 
3.899 
0.0057 
0.0049 
0.0050 


Table 4 Pairwise Fst values (above the diagonal) and the natural logarithm of geographical 


GZL 
0.0965 


6.21 
7.06 
0.92 
6.78 
7.55 
7.30 
7.28 
7.17 
7.18 
8.01 
7.74 
7.39 
7.99 
7.82 
8.06 


distance (km) (below the diagonal) between populations of Ostrinia furnacalis 


ChF 


XZh 


HSh 


JZh 


-0.0164 0.2633 0.0560 0.1029 
0.0140 0.3246 0.0891 0.1687 
0.2525 0.0238 0.0842 


6.52 
6.38 
6.52 
7.25 
7.19 
6.93 
6.90 
7.01 
7.85 
7.31 
7.28 
7.78 
7.69 
7.94 


3.02 
6. 56 
6.61 
6.46 
6.05 
0.49 
6.82 
7.70 
7.04 
7.12 
7.4 
7.43 
7.72 


0.1701 0.1107 


GG 
0.2652 
0. 1919 
0.2130 
0.2121 


6.07 


—0.0042 0.2437 


6.86 
6.71 
6.09 
6. 14 
0.32 
7.82 
7.17 
6.% 
7.31 
7.39 
7.71 


7.20 
7.08 
6.59 
6.09 
6.09 
7.98 
7.38 
6.66 
7.03 
7.38 
7.73 


0.2711 


5.13 
6.46 
0.90 
7.19 
7.33 
6.15 
7.34 
7.01 
7.32 
7.55 


MX 
0.0549 
0.0815 
0.0207 
0. 1547 


LoY SuZ 


YAh 


0. 0660 
0.0784 
0.0242 


0.0379 0.6340 
0.0282 0.623 
—0.0055 0.6095 


0.1708 0.2434 


-0.0377 -0.0371 0.0085 
-0.0195 -0.0016 0.0977 


0.2346 


6.16 
6.78 
7.09 
7.61 
0.22 
7.22 
7.01 
7.24 
7.31 


0.2318 0.1877 


一 0.0398 0.0205 


0.0075 
6.03 
6.54 
7.79 
6.78 
6.85 
7.23 
7.14 
7.31 


5.64 
7.94 
7.12 
6.36 
7.39 
7.09 
7.32 


0.4747 
0.6171 
0.6278 
0.2328 
0.6151 
0.6127 
0.5832 


8.04 
7.29 
3.80 
7.49 
7.08 
7.33 


Chj 
0.1872 
0.2523 
0.1649 
0.0778 
0.0451 
—0.0201 
0.2975 
0.0302 
0.0542 
0.1719 
0.6358 


7.64 
8.12 
7.84 
8.12 
8. 18 


PX 


DY KM 


0.0057 
0.0060 0.0048 


GZh 
2. 865 
2. 835 
3.007 
5.319 
3.165 
3.219 
0.403 
3.096 
3.151 
3.261 
10. 894 
3.406 
3.481 
3.330 
3.865 


亚洲 玉米 晶 成 对 种 群 间 的 固定 系数 (Fest 值 ) (上 三 角 ) 与 地 理 距 离 (km) 的 自然 对 数值 (下 三 角 ) 


GZh 


0.0531 0.1858 
0.0006 0.2621 
-0.0013 0. 1645 
0.1766 0.0833 
0.0040 ”0.0391 


0.0294 0.0115 
0.0083 0.01% 


0. 1668 
0.0156 


0.3182 
0.0702 


0.0450 


—0.0179 0.0658 


0. 1569 


0.1434 0.3067 
—0.0065 0.0251 
—0.0048 0.0450 
—0.0008 0. 1631 
0.5534 0.6420 
0. 0993 
0. 1042 

7.36 

6.46 

7.14 

7.32 


7.48 
6. 88 
7.41 


0. 1022 
0.0120 
0.0103 


0.2321 
0.0750 
0.0715 


—0.0174 0.0029 


0.5031 0.6278 


-0.0201 0.1206 0.2469 
—0.0276 0.0305 


0.1234 0.2468 
0.0109 
6.98 


6.94 6.% 


SYa 
3.295 
3.209 
3. 390 
5.241 
3.309 
3.383 
6.469 
3.269 
3.289 
3.457 
10. 540 
3.300 
3.711 
3.401 
4.090 
3.201 


SYa 
0.0379 
0.0282 


—0.0090 -0.0212 -0.0055 


0.2434 
0. 0085 
0.0977 
0. 1877 
0. 0205 
0.0075 
—0.0385 
0. 5832 
0. 1719 
一 0.0008 
0. 1631 
—0.0174 
0.0029 


通过 AMOVA 分 子 变 异 分 析 结 采 可 以 看 出 ( 表 ”方差 组 分 为 1.7411, 占 方差 比率 的 76.45%, 且 差 


5) ,中国 亚洲 玉米 蝇 17 个 地 理 种 群 之 间 的 方差 组 
分 为 0.5363, 占 方差 比率 的 23.55% ; 种 群 内 部 的 。” 分 化 主要 来 日 种 群 内 部 。 


异 显 车 (P< 0.05)。 说 明 中 国 亚洲 玉米 蜡 的 遗传 
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为 了 检测 亚洲 玉米 蜡 不 同 种 群 之 间 的 遗产 分 化 
程度 是 否 与 地 理 隔 离 有 关 , 通过 ZT 软件 包 进 行 
Mantel 检验 ,计算 成 对 种 群 间 的 遗传 距离 | Fst/(1 

-Fst) ] 和 矩阵 和 各 采集 地 点 间 的 地 理 距 离 (km) 的 自 


然 对 数 矩 阵 之 间 的 相关 系数 ,得 到 遗传 距离 与 地 理 
距离 之 间 的 相关 系数 > 为 -0.136 (P=0.1134> 
0.05, 10 000 次 随机 抽样 ), 表明 两 者 之 间 未 呈现 
出 显著 的 相关 性 。 


表 5 亚洲 玉米 蝇 17 个 地 理 种 群 线粒体 DNA CO 了 I 基因 序列 变异 的 分 子 变异 分 析 
TableS Molecular variance (AMOVA) analysis of mitochondrial DNA CO I gene sequences in 17 geographic 


populations of Ostrinia furnacalis 


变异 来 源 自由 度 平方 和 

Source of variation df Sum of squares 
种 群 间 Among populations 16 236. 265 
种 群 内 Within populations 396 689. 469 
总 变异 Total 412 925.734 


因 定 系数 Fixation index, Fst: 0.2355. 
入 § 
3 讨论 


一 个 物种 对 环境 变化 的 适应 能 力主 要 取决 于 其 
群体 内 部 的 遗传 多 样 性 和 相应 的 遗传 结构 。 物 种 内 
基因 流 的 模式 与 强度 ,及 其 与 种 群 遗传 结构 的 关 
系 , 是 重要 的 进化 生态 学 问题 ( 祖 元 刚 等 ，1999 ) 。 
种 群 间 的 遗传 差异 主要 是 突变 、 遗 传 漂 变 、 选 择 和 基 
因 流 等 因素 相互 作用 的 结果 。 其 中 , 突变、 遗传 漂 
变 和 选择 压力 会 促进 种 内 遗传 分 化 的 产生 ; 而 基因 
流 则 通过 配子 、 个 体 或 整个 群体 的 迁移 , 使 种 群 间 
保持 遗传 上 的 相似 性 ， 弱 化 群体 间 的 遗传 差异 
( Whitlock and McCauley，1999 ) 。 目前 对 群体 间 基 
流 的 研究 ， 主 要 是 通过 考察 不 同 群体 间 的 等 位 基 
因 频 率 或 DNA 序列 间 的 差异 来 实现 。 

在 玉米 坚 种 群 分 化 及 基因 流 方面 , 国内 外 已 有 
一 些 相 关 人 研究 报道 。 如 法 国 的 欧洲 玉米 旦 Osirinia 
nubilalis( Hiibner) 因 寄主 植物 分 化 ,致使 不 同 食性 
的 种 群 产 生 了 遗传 分 化 , 并 形成 生殖 隔离 ,两 品系 
之 间 的 基因 流 受 到 很 大 程度 的 阻碍 ( Martel ei al.， 
2003 ) 。Coates 等 (2004 ) 通过 分 析 欧 洲 玉 米 量 的 CO 

I 和 COLI 基因 序列 ， 发现 北美 欧洲 玉米 虹 不 同 地 
理 种 群 和 不 同化 性 种 群 均 存 在 着 遗传 分 化 。 在 国内 
也 有 关于 亚洲 玉米 昌 种 群 分 化 的 初步 俩 究 ， 孙 姗 等 
(2000) 用 RAPD 的 方法 研究 发 现 中 国 亚 洲 玉 米 晓 5 
个 地 理 种 群 之 间 已 经 产生 了 一 定 程度 的 遗传 分 化 。 
刘 衬 等 (2003 ) 同样 利用 RAPD 分 子 标记 技术 , 对 中 
国 不 同 玉米 种 植 区 亚洲 玉米 螟 13 个 地 理 种 群 进行 
遗传 差异 分 析 , 表明 其 种 群 间 的 遗传 差异 和 地 理 位 
置 存在 一 定 的 相关 性 , 推测 中 国 亚洲 玉米 螟 的 遗传 


方差 组 分 方差 比率 P 值 
Variance components Percentage of variation P value 
0.5363 Va 23.55 <0.0001 
1.7411 Vb 76.45 <0.0001 
2.2774 100 


分 化 主要 是 由 于 地 理 隔离 阻碍 了 基因 交流 而 产生 。 

本 研究 采用 线粒体 DNA 分 子 标记 的 方法 ,对 
中 国 亚洲 玉米 蜡 不 同 地 理 种 群 的 遗传 分 化 和 基因 流 
进行 了 更 为 系统 、 全 面 的 研究 。 据 报道 欧洲 玉米 蜡 
作为 外 来 入侵 种 , 在 20 世纪 初 才 由 欧洲 传人 美国 ， 
北美 的 欧洲 玉米 蜡 种 群 经 历 了 传人 选择 或 线粒体 
基因 的 遗传 漂 变 等 遗传 瓶颈 效应 ,因此 分 布 在 北美 
的 欧洲 玉米 量 种 群 线粒体 基因 多 态 性 水 平 很 低 
(Coates et al., 2004 ) 。 而 本 研究 结果 发 现 , 中 国 亚 
洲 玉 米 蜡 的 mtDNA CO 本 基因 存在 34 个 多 态 性 位 
点 共 35 种 单 倍 型 ,总 群体 的 单 倍 型 多 样 性 指数 Hd 
为 0.811, 种 群 内 的 单 倍 型 多 样 度 平均 为 0.712, 表 
明 中 国 亚洲 玉 米 蜡 的 线粒体 基因 具有 较 高 的 多 态 
性 , 推测 这 主要 是 由 于 亚洲 玉米 蜡 在 玉米 引入 中 国 
之 前 以 合子 和 其 他 杂 草 为 寄主 , 在 中 国 进化 的 历史 
比较 悠久 ,因而 具有 非常 丰富 的 遗传 基础 。 单 倍 型 
多 样 度 和 核 背 酸 多 样 度 在 中 国 不 同 亚洲 玉米 蜡 种 群 
中 存在 差异 , 其 原因 可 能 是 由 于 亚洲 玉米 曝 在 中 国 
的 分 布 范 围 干 分 广泛 ， 其 不 同 种 群 所 处 的 目 然 生态 
环境 不 同 , 在 自然 选择 及 寄主 植物 的 选择 作用 下 ， 
种 群 的 遗传 结构 已 发 生 了 较 大 变异 。 

固定 系数 Fst (或 称 基 因 分 化 系数 ) 作为 遗传 
距离 的 一 种 常用 测度 参数 , 可 反映 出 不 同 物种 或 种 
群 之 间 的 基因 差异 程度 。 中 国 亚洲 玉米 螟 不 同 地 理 
种 群 之 间 的 Fsi 值 在 -0.0398 ~ 0. 6420 之 间 , 总 群 
体 的 Fst 值 为 0.234。 据 报道 , 基因 分 化 系数 的 值 
介 于 0 ~0.05 之 间 表 示 群 体 间 遗 传 分 化 程度 很 弱 ; 
0.05 ~0.15 为 分 化 程度 中 等 ; 0. 15 ~ 0. 25 为 分 化 
程度 较 大 ; 大 于 0. 25 为 分 化 程度 很 大 ( Rousset， 
1997)。 因 此 从 总 体 上 来 看 , 不 同 地 理 种 群 之 间 的 
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遗传 分 化 程度 较 高 。 通 过 中 性 检验 分 析 结 有 果 可 以 推 
测 ， 中国 的 亚洲 玉米 坚 种 群 在 过 去 没有 出 现 群 体 扩 
张 和 持 续 增 长 模式 , 群体 大 小 保持 相对 稳定 状态 。 

基因 流 是 种 群 遗 传 结构 均 质 化 的 主要 因素 之 
一 ， 具 有 高 水 平 基因 流 的 物种 往往 比 具 有 有 限 基 因 
流 的 物种 遗传 分 化 小 ，Whitlock 和 McCauley(1999 ) 
认为 当 Nm >4 时 , 种 群 间 的 基因 交流 比较 充分 , 匀 
质 化 作用 足以 抵制 遗传 漂 变 的 作用 , 防止 种 群 间 遗 
传 分 化 的 产生 。 本 研究 中 , 由 Nei 氏 指 数 计算 的 总 
群体 基因 流 Nm =2. 02 <4, 表明 种 群 间 的 基因 流 较 
小 , 这 与 玉米 旦 群体 间 遗 传 分 化 程度 较 高 的 结论 相 
对 应 。 同 时 也 反映 出 ,虽然 亚洲 玉米 虹 成 虫 扩散 能 
力 较 强 ( 王 振 营 等 , 1994), 但 它们 在 不 同 地 理 种 群 
间 的 迁移 扩散 机 会 相对 有 限 ， 从 而 限制 了 种 群 之 间 
的 基因 交流 。 这 一 结果 为 今后 制定 转基因 抗 虫 玉米 
的 抗 性 治理 策略 中 庇护 所 的 设置 范围 和 距离 提供 了 
科学 依据 。 

通过 对 单 倍 型 在 不 同 种 群 中 的 分 布 情况 及 各 单 
倍 型 的 系统 进化 分 析 , 结果 发 现 单 售 型 Hap3 是 17 
个 亚洲 玉米 螟 种 群 的 共享 主体 单 倍 型 ,其 可 能 是 玉 
米 蝶 进化 历史 中 出 现 的 较为 稳定 、 环 境 适应 性 较 强 
的 一 种 单 倍 型 。 从 各 单 倍 型 在 邻接 树 (NJ 树 ) 中 的 
分 布 和 各 单 倍 型 在 实际 地 理 种 群 中 的 分 布 情况 来 
看 , 单 倍 型 之 间 和 其 分 布 的 地 理 单 元 之 间 并 没有 明 
显 的 对 应 关系 , 各 地 理 种 群 具 有 的 单 倍 型 都 处 于 相 
互 散布 的 、 比 较 混 杂 的 分 布 格局 中 。 男 外 ,通过 各 
种 群 间 和 遗传 距离 及 地 理 距离 的 相关 性 分 析 ，, 未 发 现 
二 者 具有 显著 的 相关 性 。 因 此 , 地 理 隔离 并 不 是 造 
成 中 国 亚洲 玉米 蜡 种 群 间 遗传 分 化 的 主导 因素 。 关 
于 亚洲 玉米 螟 在 中 国 的 地 理 遗 传 分 布 格局 的 演化 形 
成 的 原因 ， 有 待 于 今后 进一步 深入 研究 。 
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